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A quasi-analytical model for energy-delay-reliability 
tradeoff studies during write operations in 
perpendicular STT-RAM cell 

Kamaram Munira, William H. Butler and Avik W. Ghosh 



Abstract — One of the biggest challenges the current STT-RAM 
industry faces is maintaining a high thermal stability while 
trying to switch within a given voltage pulse and energy cost. 
In this paper, we present a physics based analytical model that 
uses a modified Simmons' tunneling expression to capture the 
spin dependent tunneling in a magnetic tunnel junction(MTJ). 
Coupled with an analytical derivation of the critical switching 
current based on the Landau-Lifshitz-Gilbert equation, and the 
write error rate derived from a solution to the Fokker-Planck 
equation, this model provides us a quick estimate of the energy- 
delay-reliability tradeoffs in perpendicular STTRAMs due to 
thermal fluctuations. In other words, the model provides a simple 
way to calculate the energy consumed during a write operation 
that ensures a certain error rate and delay time, while being 
numerically far less intensive than a full-fledged stochastic cal- 
culation. We calculate the worst case energy consumption during 
anti-parallel (AP) to parallel (P) and P to AP switchings and 
quantify how increasing the anisotropy field Hk and lowering the 
saturation magnetization Ms, can significantly reduce the energy 
consumption. A case study on how manufacturing variations of 
the MTJ cell can affect the energy consumption and delay is also 
reported. 

Index Terms — Spin polarized transport, Spin torque 
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Fig. 1. Three coupled equation blocks allow us to compute the energy 
consumption for a write operation in a perpendicular STTRAM cell within a 
given switching time and error threshold. Equation block 1: Given a write- 
error-rate (WER) and delay time r, we use the non-switching expression, Eq.l, 
to extract the overdrive I/Iq ratio, where I is the current applied and Iq is 
the critical switching current. Equation block 2: gives us the critical switching 
current Iq using Eq. 2. From these two blocks, we get /. Equation block 3: 
Using a modified version of Simmons equation, Eq. 3, we can then estimate 
the switching voltage V and thence the switching energy IV r. The last two 
steps have to be solved self-consistently, as the switching voltage for a given 
WER depends on current Iq , which in turn depends on the voltage-dependent 
polarization 77 (V), which further depends on the transport parameters Ep, 
m c and 5. 



I. Introduction 

WITH its fast write and read, small cell size, non- 
volatility and excellent endurance, Spin Transfer 
Torque-based Random Access Memory (STT-RAM) devices 
bear excellent potential of dominating the embedded and 
standalone memory world in the near future Q, Q. However, 
in order to compete with existing embedded technology, an 
STT-RAM bit cell should be able to meet the following 
requirements: (i) a tunneling magnetoresistance (TMR) needs 
to exceed 150% for low-power read; (ii) thermal stability (A) 
must be greater than 60 to achieve high static reliability; (iii) 
a scalable spin transfer torque switching current (J c < 1.5 x 
10 6 A/cm 2 ); (iv) a very low write error rate (WER) ~ 10 -9 
with a fast device write time < 10ns; and (v) an overall low 
energy write (<lpJ), which corresponds to a low switching 
voltage 1 3 ]. Achieving such a diverse set of targets will require 
meticulous design and material engineering. 

The reliability of the writing process in arrays of STT-RAM 
cells is dependent on the following factors: (i) variation of 
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the electronic, magnetic and geometrical parameters across the 
array (Fig. 1)0, Q, (ii) thermally activated initial angle, 6q, at 
the beginning of the writing process, (iii) thermal fluctuations 
during the writing process! 6] and (iv) voltage drop across the 
1 -transistor 1-MTJ cell[7|. For this paper, we disregard the 
voltage drop across the transistor in (iv) as it it outside the 
scope of this work. We address the rest in our model. 

(ii) and (iii) are the dominant players when it comes to 
the writing process in STT-RAM, with (ii) being the most 
critical. The free layer in the STT-RAM MTJ memory cell is 
responsible for data storage, and must have enough stability 
against stochastic thermal switching to retain data for at least 
10 years Q. In other words, the thermal stability, A, needs 
to be greater than 60. While we desire the device to have 
high thermal stability, the critical current and the time needed 
to switch also increase proportionally. Higher thermal stability 
also contributes to longer switching tails, as the magnetization 
tends to be nearer to the stagnation points (0 = or tt) 
in the potential landscape, without any external perturbation. 
Near the stagnation point, the torque proportional to sin# is 
too small to initiate switching. Therefore a thermally stable 
free layer would require a relatively high current to have very 
low error rate 10 ~ 9 for switching within 5-10 ns, resulting in 
costlier write operations. 

While full-fledged stochastic material and circuit Simula- 
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tions can help narrow down the target compositions, they 
are time-consuming. Therefore, it is important to develop in 
parallel a fast, quasi-analytical physics based model to explore 
the multi-parameters optimization problem involved. This is 
the aim of the paper. 

In this paper, we introduce a quasi-analytical model (Fig 1) 
that can be used to calculate the current and energy consumed 
during a write operation with an acceptable WER for a 
given pulse width, r, in a perpendicular STT-RAM cell. The 
calculation is done in three blocks: (1) Given a write-error- 
rate and switching time, we use an expression for switching 
errors from [ 8 ] to extract the I/Ic ratio, where / is the current 
applied and Ic is the critical current. (2) The critical switching 
current Ic is calculated using the Landau-Lifschitz-Gilbert 
(LLG) equation to yield the applied current /. This expression 
also includes additional effects due to voltage-dependent po- 
larization responsible for observed switching asymmetries |9l. 
(3) Using a modified version of Simmons tunneling current 
equation, thats includes the crucial pre-factors responsible for 
TMR, we can then estimate the switching voltage and thence 
the switching energy cost. The last two steps are coupled and 
need to be solved self-consistently, as the switching voltage 
depends on critical current, which depends on polarization, 
which in turn depends back on voltage. The energy consumed 
during the write operation can then be estimated from the 
switching current, voltage and time as IV r. 

Coupling all three equations - the WER expression, the 
modified Simmons tunneling equation and the modified critical 
current equation! 10] including the role of voltage-dependent 
polarization are non-standard in the literature. Note that this 
paper focuses on perpendicular materials, for which the WER 
expression has been explicitly worked out and calibrated 
against stochastic LLG solutions! 8]. In future papers, we 
will extend this study to in-plane materials as well as to 
double barrier STT-RAMs. Section II describes the worst case 
current needed for perpendicular STT-RAM, based on solving 
the Fokker-Planck equation. In section III, we present the 
Simmons' equation modified to include spin dependent tun- 
neling. The modified equations are used to extract electronic 
parameters from published CoFeB/MgO/CoFeB MTJ current- 
voltage graphs. Using the parameters extracted from the fit in 
section III, the worst case energy consumption is thereafter 
calculated for a given A and error rate in section IV. The 
effect of manufacturing variations on energy consumption is 
also studied. 




10 15 
Delay, x (ns) 



Fig. 2. Write error rate vs. delay time at room temperature for A =60. To 
attain an error rate of 10 -9 within given delay times 5ns, 10ns and 20ns , 
current pulses of 3.232/c, 2.104/c and 1.543/c are needed. We use Hk 
of 3.34kOe and Ms of 1257.3 emu/cc for this estimate fTTl . 



the simple case of a perpendicular material with uniaxial 
anisotropy in the macrospin approximation. The corresponding 
probability of not switching within a given delay r, defining 
the WER, is expressed as 

-7T 2 A(i - l)/4 



WER(t) = 1 - exp 



a2a 1 H K T{i-l)/{l+a' 2 ) _ l 



c 



(1) 



where / is the current required and Ic is the critical switching 
current, a is the Gilbert damping parameter, 7 is the gy- 
romagnetic ratio, and A = HKMs^t/2kBT is the thermal 
stability of the free layer, ft is the volume of the free layer 
and Ms is the saturation magnetization. Hk = H^ — AirMg 
is the effective anisotropy field resulting from both crystalline 
anisotropy (H^) and demagnetization field. HkMs^I/2 is the 
energy barrier that separates the two magnetization directions, 

= and tt. Eqn. 1 is valid for precessional switching where 

1 » 1. 

Solving the LLG equation we get a critical current similar 
to the expression derived by J.Z. Sun fTOi but modified to 
include the proper switching asymmetry through the effective 
voltage dependent polarization(described in the next section), 



I c = 2aeH K M s n/hri(V) 



(2) 



From the given WER and r, we can get the current overdrive 
i = I/Ic- In conjunction with the critical current above, 
we can then extract the applied current /. Finally, using the 
modified version of Simmons equation described in the next 
section, we can then estimate the switching voltage. 



II. Current needed for worst case writing due to 

THERMAL FLUCTUATIONS 

Given a desired write-error-rate (WER) and switching delay 
r, we can calculate the amount of current overdrive required 
relative to the critical switching current 0. Numerically this 
amounts to introducing a stochastic Langevin thermal torque 
in LLG and then taking a time average over many runs. The 
equivalent probabilistic approach is to solve the Fokker-Planck 
equation for the probability density of each configuration and 
then average over the distribution to extract various moments. 
The Fokker-Planck equation can be solved analytically for 



III. Modified Simmons' equation for magnetic 

TUNNEL JUNCTION 

In a magnetic tunnel junction, there are two ferromagnetic 
(FM) electrodes with an insulator between them. Ignoring 
atomistic effects such as symmetry filtering! 12], the one-band 
alignment problem across the magnetic tunnel junction can 
be specified with six parameters: (i) U is the barrier offset 
between the contact and the insulator, (ii) Ep is Fermi-energy 
relative to the conduction band minimum in each contact 
(assumed to be comprised of the same material), (iii) 5 is the 
band splitting between majority and minority spin electrons 
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Fig. 3. Schematic band alignment diagram of an MTJ. The bottoms of the 
t and J, majority conduction electron bands in the FM contacts are labeled 
E™ aj and E™ in respectively, and are separated by 5 eV, while the insulating 
oxide introduces a tunnel barrier. Ep is the Fermi-energy relative to the 
conduction band minimum, W is the width of the insulating barrier and U is 
is the barrier offset between the contact and the insulator. 



Fig. 4. Parameter extraction from published CoFeB/MgO/CoFeB MTJ 
resistance- voltage curve 1 16 1 with modified Simmons' model. The parameters 
extracted from the model are Ep =2.2eV, U= 1 eV, W = 1 nm, m c = 0.3 
rao> mbarr =0.18 mo and 8 = 1.98 eV. Curve A is for anti-parallel mode 
while D for parallel mode, mo is the mass of an electron. 



and related to the contact polarization, (iv) W is the width of 
the insulator, while (v) m c and (vi) mfo arr are the effective 
electron masses in the contact and barrier respectively (Fig. 
3). At thermodynamic equilibrium, the Fermi levels at the 
two FM electrodes are aligned. However for a positive bias 
qV applied to the free layer, the energy levels in the fixed 
FM are shifted up by qV/2 and the ones in the free FM 
are lowered by qV/2, separating the Fermi levels at the two 
electrodes by qV. The mismatch in the Fermi levels cause 
current to flow in the system through oxide tunneling. In 
a free electron model for the electrodes, the longitudinal 
spin-polarized electron momentum at a given energy E in 
each of the three regions can be written as follows. Fixed 
FM: k tA (V) = y / (2m c [E F - (5 T 6 - qV)/2])/h, Insulator: 
k = \J 2vrt\ )arr U I ft where (3 — KVfi c j vrt\) arr , and Free FM: 
k+^V) = ^(2m c [E F — (5 =F 5 + qV)/2])/h. The magneti- 
zation of the left contact is fixed, while the magnetization of 
the right contact is free to rotate and is defined by the relative 
angle, 6. 

Using the methodology used in [13], a weighted sum 
of the tunneling up and down spin electrons was obtained, 
with weighting factors given by their respective densities 
of states. The transverse modes with momentum fey, and a 
linear potential drop in the insulator were included using the 
WKB approximation prescribed by Simmons in [14]. The total 
current density, J, at a given voltage and angle 6 in a magnetic 
tunnel junction is, 



Alt 2 TlW 2 



2W 

(U -qV/2)e~~h~ 
2W 

-(U + qV/2)e h~ 
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From the spin current, we can calculate the torque 
exerted on the free magnetic layer fT5lL T(V) = 
(h/2q)J avg (V)r}(V)[m x (m x m p )] where J avg (V) = 
[J(V,6 = 0) + J(V,0 = tt)]/2. m and rh p represent the 
unit vector orientations for the magnetizations in the free 
and fixed layers respectively. When a positive bias is applied, 
polarized electrons flow from the fixed FM contact to the free 
layer, causing the magnetization of the free layer to switch 
from anti-parallel to parallel. For negative bias, the electron 
flow removes magnetization from the right layer, causing it 
to rotate from P to AP. Notice that the torque exerted at 
the right, free ferromagnetic layer depends on the effective 
voltage -dependent polarization of the incident electrons from 
the left, fixed layer[9|. This voltage dependence is critical to 
understanding the observed P to AP vs AP to P asymmetry in 
STTRAMs. The effective polarization of the incident electrons 
at a given voltage is given by, 



v(v) 



kf 



fcf- + ki 



(4) 



J(V,0) = [J t (V,9) + Jj r (V,O)]J o (V) 



(3) 



The voltage dependency on electron momentums, k a and k+ 
are not explicitly mentioned in Eq. 3 and 4 for notational 
simplicity. Fig. 4 show results from the analytical model 
fitted with a published experiment. Eq. 3 was used to fit the 
resistance vs voltage data for in-plane CoFeB/MgO/CoFeB 
MTJ fToll . The parameters extracted from the model are listed 
in the figure caption. 

IV. Worst case energy consumption and switching 

DELAY DUE TO THERMAL FLUCTUATIONS AND 
MANUFACTURING VARIATIONS 

Using the self-consistent analytical model illustrated in Fig. 
1, we calculate the worst case energy consumption during a 
write operation at room temperature for A=60 and various de- 
lays (Fig. 5) due to thermal fluctuations. Energy consumption 
during P to AP switching is greater than that consumed during 
AP to P switching. This is because the effective polarization 
of the torque at negative voltage is less than that at positive 
voltage (9|, (331 • However, for a magnetic tunnel junction with 
half-metallic contacts (<5=3eV), energy consumption during AP 
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Fig. 5. Worst case energy consumption during a write operation at room 
temperature for A=60 and error rate of 10 -9 . Energy consumption during (a) 
AP to P switching is less than that consumed during (b) P to AP switching. 
This is because the effective polarization of the torque at negative voltage is 
less than that at positive voltage. Energy consumed is the product of current, 
voltage applied and delay (IV t). At lower delay, the current, and hence 
energy, decreases logarithmically with delay. At higher delay, the energy 
consumption increases linearly with delay. 



to P switching is equal to that during P to AP switching 
because the effective polarization of the torque at negative 
voltage is same as that at positive voltage (Fig. 6). 

Energy consumed is the product of current, voltage applied 
and delay (IV r). The voltage drop considered here is the one 
across the MTJ. The additional voltage drop in the transistor 
which usually accompanies the STT-RAM cell is disregarded 
in this paper. When the applied current is greater than the 
critical current, the switching is precessional and the switching 
delay of the free layer is inversely proportional to the current, 
I oc I c + Ic(X + a 2 )ln(7r/20 )/(ajH K r) cx (At + ff)lfl8l 
[8], where #o is the initial angle of the switching free ferromag- 
net. Assuming that the current in the MTJ grows exponentially 
with applied voltage (Eq. 3), V cx [ln(I / C) / D) 2 . The energy 
consumption with respect to delay is 



Energy cx IV t cx (At + B) 



^ Ct) 



D 



E (5) 



where A(=J C ), B(=/ c (l+a 2 )/n(7r/2<9 ) /(ajH K )), C(=Area- 
q 2 /4ir 2 hW 2 ), V(=2W^/2jn~ c /h) and E are constants. The 
units of constants A, B, C and E are Amp, Amp- sec, Amp, 
V~ x l 2 and J respectively. At lower delay, the energy de- 
creases logarithmically with delay as [ln(I/C)/D] 2 is the 
dominant term. At higher delay, the energy consumption 
increases linearly with delay as At is the dominant term. 
There is thus an energy minimum, E m i n , corresponding to an 
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Fig. 6. Keeping volume and thermal stability of the free layer the same 
(A=60), the minimum energy consumption can be reduced by increased 
the anisotropy field, Hk and decreasing the saturation magnetization, Ms 
accordingly. The duration of the voltage pulse needed to achieve the switching 
with 10 _ 9 error rate also decreases with higher anisotropy field a=0.02. 
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Fig. 7. Worst case energy consumption during a write operation at room 
temperature for A=60 and error rate of 10 _ 9 for a magnetic tunnel junction 
with half-metallic contacts (5=3eV). Energy consumption during AP to P 
switching is equal to that during P to AP switching because the effective 
polarization of the torque at negative voltage is same as that at positive voltage. 
H K = 5.34/cOe and M s = 7 86. 5 emit /cc. 



optimal delay value, T m i n . Fig. 5 shows a fit of the numerical 
results with Eq. 5. The constants to fit the energy-delay 
curve for AP to P switching at Hk = 3.34/cOe were: 
A=7.2mA, B=1.8xlO- n A-s, C=5.7mA, D=3.83V~ 1 / 2 and 
D=0.5xl(T 14 J. 

In addition, keeping A at 60, if Hk is increased from 
3.34 kOe to 5.34 kOe, and Ms adjusted accordingly from 
1257.3 emu/cc to 786.5 emu/cc, the E m i n is shifted from 
0.1252pJ at 10.5ns to 0.0783pJ at 6.6ns for AP to P switching. 
Similarly for P to AP switching, the energy minimum shifts 
lower from 0.3829pJ at 16.25ns to 0.2393pJ at 10ns. In 
other words, the E min can be reduced quite significantly 
by increasing the anisotropy field, Hk and decreasing the 
saturation magnetization, Ms accordingly, while maintaining 
thermal stability. The duration of the voltage pulse needed to 
achieve the switching with 10 ~ 9 write error rate, r m i n , also 
decreases with higher anisotropy field (Fig. 7). 

All the above simulation results (Figs. 5-7) only considered 
thermally activated initial angles, #o, at the beginning of the 
writing process and thermal fluctuations during the writing 
process in a given MTJ cell with a fixed dimension, electronic 
and magnetic parameters. However, individual MTJ cells in 
STT-RAM arrays will not all have the same physical dimen- 
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Fig. 8. Effect of manufacturing variations on minimum energy consump- 
tion, E m i n , during a P to AP write operation with an error threshold of 
10 -9 . The varying parameters are effective anisotropy field (Hk), saturation 
magnetization (Ms), Gilbert damping parameter (a), cross-sectional area of 
the MTJ, thickness of the free layer (t), width of the tunneling barrier (W), 
and ferromagnet-oxide interface (U). While variation in the thickness of the 
tunneling barrier has the most effect on the delay, cross-sectional area has 
the least. Mean value of the MTJ parameters used are Hk = 5.340e, 
Ms = 785.6emu/cc, a = 0.02, radius of the free layer=20nm, t=2nm, 
W=lnm and U=leV. 

sion or parameters because of manufacturing variations. The 
key parameters, Hk, Ms, ol, cross-sectional area of the MTJ, 
thickness of the free layer (t), W and U can vary by 10% 
from their mean value |4'|. The self-consistent analytical model 
illustrated in Fig. 1 can be used to evaluate the write perfor- 
mance against various intrinsic variabilities as show in Figs. 8 
and 9. All the above parameters were varied individually by 
±10% and their effect on the E min and r min are reported. 
The error rate is fixed at 10 -9 . Differing Hk, Ms, t and 
cross-sectional area will affect A. Increased A and a increases 
the critical current and energy consumption. Increasing barrier 
width and height results in more voltage being needed to get 
the desired current. While variation in the thickness of the 
tunneling barrier has the most effect, cross-sectional area has 
the least. The worst-case energy consumption and delay in an 
STT-RAM array can be calculated by considering an MTJ cell 
with biggest variation (+10%) in all the parameters above. The 
values are 1.03pJ at 5.25ns for P to AP writing. A detailed 
study of the manufacturing variations will be reported in a 
later paper. 

V. Conclusion 

In this paper, we present two models: (1) a modified 
Simmons' tunneling expression that precisely captures the spin 
dependent tunneling in a magnetic tunnel junction (Eq. 3) 
and (2) a corresponding model that approximately calculates 
the energy consumed during a write operation that ensures 
a certain error rate and delay time. The combination of the 
two models enables us to calculate the energy consumption 
without performing a full-fledged stochastic calculation and 
to evaluate ways to reduce the energy consumption during 
worst-case writing in STT-RAM. At lower delay, the energy 
consumption decreases logarithmically with delay while at 
higher delay, it increases linearly with delay, creating a delay 
point where the energy consumption is minimum (Fig. 5). We 
also quantified how increasing the anisotropy field Hk and 
lowering the saturation magnetization Ms while keeping A 
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Fig. 9. Corresponding minimum switching delay during a P to AP write 
operation with an error threshold of 10 -9 . The varying parameters are 
effective anisotropy field (Hk), saturation magnetization (Ms), Gilbert 
damping parameter (a), cross-sectional area of the MTJ, thickness of the 
free layer (t), width of the tunneling barrier (W), and ferromagnet-oxide 
interface(U). 

fixed, can significantly reduce the energy consumption and 
the switching delay. The effect of the variations of geometric, 
magnetic and electronic properties in STT-RAM arrays is 
studied and it is seen that variation in the thickness of the 
oxide has the biggest role in determining the worst case energy 
consumption scenario. 
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